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Abstract 

It is shown that when the Einstein-Hilbert Lagrangian is considered without any non-covariant 
modifications or change of variables, its Hamiltonian formulation leads to results consistent with 
principles of General Relativity. The first-class constraints of such a Hamiltonian formulation, with 
the metric tensor taken as a canonical variable, allow one to derive the generator of gauge trans- 
formations, which directly leads to diffeomorphism invariance. The given Hamiltonian formulation 
preserves general covariance of the transformations derivable from it. This characteristic should 
be used as the crucial consistency requirement that must be met by any Hamiltonian formulation 
of General Relativity. 
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I. INTRODUCTION 



The Hamiltonian formulation of General Relativity (GR) is more than a half century old. 
To cast GR into Hamiltonian form, a method which deals with singular systems is required. 
A way of dealing with this problem, the constrained dynamics, was proposed by Dirac in 
1949 [1]. Dirac's procedure was almost immediately applied by Pirani, Schild, and Skinner 



to GR using the covariant metric tensor as the canonical variable. However, the analysis 
of [2I was not complete; the time development of secondary constraints was not considered 
and, strictly speaking, the closure of the Dirac procedure was not demonstrated. Neither the 
Dirac conjecture that all first-class constraints generate gauge symmetries [3j] nor explicit 
methods for derivation of gauge transformations were known at that time. 

A few years later, Dirac revisited this problem himself J4] with the same choice of a 
canonical variable, the covariant metric tensor, and made a modification of the Einstein- 
Hilbert (EH) Lagrangian to simplify the primary constraints. This modification does not 
affect equations of motion but, according to Dirac, "can be achieved only at the expense of 
abandoning four- dimensional symmetry" [4]. He explicitly demonstrated that the canonical 
Hamiltonian is proportional to a linear combination of secondary constraints. However, as 
in [2I, it was not possible at that time to derive the gauge invariance which results from the 
presence of first-class constraints. 

Shortly after the appearance of Dirac's article [4|, a new set of variables was introduced 
by Arnowitt, Deser, and Misner (ADM) (see [s| and references therein) and the geometrical 



interpretation of these new variables was widely emphasized and discussec. 
The Dirac conjecture was converted into an algorithm 3] (see also 



e.g. see 
9j) long after the 



appearance of . It was applied only to the most widely accepted formulation based 

on the ADM variables. The complete derivation of the gauge transformations from the full 

n 

set of first-class constraints in the ADM formulation actually appeared only recently [10J. 



However, a field- dependent redefinition of the gauge parameters found in 
to present the result in the form of diffeomorphism. These transformations 



known for a long time 



Ofl is required 



1 11 ] and were partially derived in [7J for the ADM formulation as 



101 ] have been 



an illustration of a general procedure for the derivation of gauge transformations. In 11] 
the distinction between these transformations and the diffeomorphism transformation was 
pointed out. In [10|] this distinction was called "the unity of the different symmetries" (our 
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Italic). 

Why is it in the case of GR, where the Lagrangian and the equations of motion are 
invariant under the diffeomorphism transformation 1 

8(diff)9nv = —£,ti;v — £,v;n , (1) 

does the Hamiltonian formulation lead to "unity" , instead of giving the equivalent result? 
The Lagrangian and Hamiltonian formulations are supposed to give the same results; in 
particular, the same gauge invariance. This is well-known equivalence for ordinary gauge 



field theories. In addition, we refer to the result of Samanta 15j where ([T]) was derived using 



16| without introducing any new variables or modifying the EH 



the Lagrangian method 
Lagrangian. 

The ADM change of variables, or rather the geometrical meaning related to it, was 
criticized by Hawking on general grounds as contradictory to "the whole spirit of relativity" 
~y\ . In our view, the disagreement between the results of 10 1 and [3| is a confirmation 
of Hawking's criticism. To eliminate such a discrepancy, we reconsider the Hamiltonian 
formulation of GR by neither using new variables (as was done by ADM) nor by making 
additional modifications of the original action (as was done by Dirac). We revert to the first 
canonical treatment of GR in [2] and prove by explicit calculations that the Hamiltonian 
formulation of GR with the metric tensor as the canonical variable leads to the symmetry 
of (CE]), as derived in the Lagrangian approach of [15|]. This is expected from a consistent 
Hamiltonian formulation of GR. 



1 In mathematical literature the term diffeomorphism refers to a mapping from one manifold to another, 
which is differentiable, one-to-one, onto, with differentiable inverse. However, in the literature on GR, the 
word "diffeomorphism" is often used as equivalent of the transformation (TT]) (the semicolon ";" means a 
covariant derivative) and in our article the latter meaning is employed. By a gauge invariance one usually 



understands the invariance in the same coordinate frame of reference |12l Il3j and with this respect the 
transformation (fT]) serves as a gauge invariance of GR and can be written in variety of forms (e.g. of 
Section V is more suitable for comparison with the results there). Note that the right-hand side of (TT]) is 
in fact a Lie derivative of the metric tensor [14j which is obviously a true tensor, so it immediately follows 
that it is not affected by any change of coordinates (in other words, (TT]) is generally covariant). 
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II. THE HAMILTON! AN 

The starting point of the Hamiltonian formulation in [2|, |4( is the r — T part 

HQ of 

the Einstein-Hilbert Lagrangian which is quadratic in first order derivatives of the metric 
tensor 



1 

4 

where 



l = V^g aP (rjy^ - r^rjjj = 7 v / =^ Q/37 ^^, 7 <w (2) 



ga/3-ypup _ yp fj,u _ a/j, f3v _|_ ^ g " P g P " g 7 _ 2g a/3 g" /fl g lyp 

To find momenta 7r Q/3 conjugate to g Q/ 3, we rewrite (T5]) to explicitly separate terms with 
time derivatives of the covariant metric tensor (the "velocities") 



L = \v=9B a ^°g af3j0 g^ + \^r g B^^ g a ^g^ k + l^gB^^g^g^ (3) 

where Latin alphabet is used for spatial components, for temporal one, and () or (... | ...) 
-brackets indicate symmetrization in two indices or two groups of indices, i.e., 

1 1 

g(oc/3)j(ivp _ _ ^a/3-ypup _|_ j^f3a~fpvp^ g(af3~/\pv 'p) _ _ ^a/3-ypup _|_ j^pupa/3-y^ 

Using (j3J) we obtain 

5L 1 , „„^ nl ..,. n \ 1 



7T 



7CT = = -'^B^'W^g^o + -V=gB^°^g^ k . (4) 



<W,o 2 _ 2 

The explicit form of the first term is (see also [4j) 



l^fl(W»M ffw0 = ^g^ET^g^ ( 5 ) 



where 



0/x Ov 

^^^-eV, e ^ = ^_y_^. ( 6 ) 

Note that both e^ u and E^ 7CT are zero, unless all //, ^, 7 and cr differ from 0, and so 
from (jSJ) it immediately follows that we have d primary constraints (d is the dimension of 
spacetime) 



-gB 



((Oo-)0|/Xi/fc) 



(7) 



For 70" = pq, equation (j4]) is invertible and giving 



9mn,0 



where l mnpq is inverse of E pqkl : 



Lmn 'Pi g oo 



(8) 



1 



z9mn9pq 9mp9nqi -Lmnpq 



'■ /t; n net 1 -' m 71 ' 



(9) 



ra " M d-2 : 

The appearance of the singularity in (J9]) for d = 2 corresponds to the fact, that in two 
dimensions Eq. (j4]) cannot be solved for all velocities. So the number of primary constraints 
equals the number of independent components of the metric tensor. 

The Hamiltonian is defined to be H = n a ^g a f3,o — L, which after elimination of the 
velocities using ([8]) gives the total Hamiltonian 



H T = H c + g 0afi ^ 
where the canonical Hamiltonian H r is 



(10) 



ff — ^ T „mn pq \_ j ran TD{pqO\auk) 

1J c / — —^QQ A mnpq'' " ^Q ± mnpq li ^ y/j,u,k 



1 



r r>((mn)0\fiuk) r)(pq0\af3t) -r>livkafit 

qq J- mnpq D D D 



9^u,kg a f3,t- 



(11) 



III. CLOSURE OF THE DIRAC PROCEDURE 



The fundamental PB are defined by 



uu 



and, following the standard procedure (see 3,Q,Q]), we have to calculate time development 
of primary constraints, {(p 0a , Ht}- Using ([7]) we now find the PB among primary constraints 



Oct ,0 Oq 



(12) 
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while, because H c is independent of 7i 0a (see (TTTj) ). {<p 0a , H c } cannot be proportional to 
primary constraints. This leads to the secondary constraints: 



1 1 n 0a 
Y 0(T = — - - - I 



.mn pq 



-± / 7r""V 

2 ^/=g g 00 npq 



(13) 



+ l^Imnpqn mn A^°^g^ k + P m Ttf + ~ \ 



K pq e km g mk 



9fj,v,k9ot/3,t 



1 / o 0<j 

I Z J j^((mn)0\fiuk) j^(pq0\a/3t) _ Oct g[iukaf3t 

- <T (2g m g q(X g tp - g m g al3 g qt + rtV - 2g«*g°Pg ot - 2g ta g 0f3 g 0q + 2g qt g 0a g^) 
+g 0a (2g f3t g ap g 0q - - 2g pq g v g 0a + 2g pt g qP g 0a + g pq g a(3 g ot - g tp g af3 g 0q )] 

- \v^gg^k9a^ [g at {g afl g Pu g ok + g^g at g op - 2<fW /3 ) 

+ g 0a {2g at g^g vk - 3g tf *g uk g a(3 - 2^VV + 0"W" + 2g^g ul3 g ka ) 
+9°" {{g aP g vt - 2g va g tp ) g ok + 2 (/V* - /VI q 0a 



where 



!T ((<f V< - W a 9 tl3 ) g ok + 2 {<P<P - g" k g tv ) g 0a + (2g kf3 g at - g^g kt ) g°»)] 



We observe that the primary constraints (|7| arc the same as those derived in [2] (sec Eq. 
(6) of S]) and provide the correct limit of linearized GR j2o|, but our (1111) is different from 
H c of B| (see Eq. (8) of jsj). While Eq. (1TT1) gives the correct limit of linearized GR, the 
H c of |2j does not. 

Of course, Eq. fflBl can be presented in different forms, however this form is convenient 
when performing an additional consistency check. Having H c proportional to the secondary 
constraints is a common feature of all generally covariant theories (e.g., see 2l|) and this 
can be seen as a good verification of the correctness of (1131) . Comparing (Hip and (1131) we 
obtain 



H c = -2g 0(TX 0a + (2go m n mk - V^9 (g 0n B^ ma ^ + g^B 



((0 7 )fc[a/9t) 



) 9aP,t) 



(15) 
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Let us continue to apply the Dirac procedure and consider the time derivatives of sec- 
ondary constraints, {x 0a , Ht}- We find that 



{x°°,r} = - 1 2 -g^x 00 



(16) 



and 



9 n a 1 1 

\X j-TCcj- — -=l m npqK ~oo"X + of #00,fcX + <>0 X,fc (,1'J 

V 9 9 £ 

# ^oo.fc + 2g na 9on,k + fi' nCT — 55- (sWi,fc + fl'fc 

This completes the proof of closure of the Dirac procedure, i.e., no further constraints 
arise. According to (112jl . (115j) . ({TBI , and ( TT71) . all constraints are first-class. This is sufficient 
to derive the gauge transformations. A number of algorithms of such a derivation are 
available s|, Q]. We follow the approach that was applied for the first time to field theories 



by Castellani 



a. 



IV. THE GENERATORS 

The Castellani procedure [?| is based on a derivation of generators of gauge transfor- 
mations which are defined by chains of first-class constraints. One starts with primary 
first-class constraint (s), i = 1,2, and from them constructs chain(s) £i Gly Here £ ( . is 
a time derivative of the zth gauge parameter of order n (n = 0, 1, ...) and n corresponds to 
the length of chain. The number of gauge parameters is equal to the number of first-class 
primary constraints. Note, that the chains are unambiguous when the primary constraints 
are defined. The functions G l ,s are calculated as follows: 

Gfo (x) = cfP* (x) , (18) 

G^ (x) = -{4P a (*) ,H T } + J a; (x, y) ^ (y) d 3 y, (19) 

where functions a° (x, y) have to be chosen so that the chain ends on the surface of primary 
constraints 
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{G? \, Ht] = primary 
and the generator G (£ CT ) is given by 



(20) 



G (£o0 — ^o-G^o) + CafiG^y (21) 

To construct the generator fl2T|) . we have to find functions a" (x, y) using the condition fl20|) 
{^ 0) ,tf T } = -{ X 0ff (x),tf T } + y {a;,// T }0°M2/)rf 3 ?/ + J a;{^(y),H T }d 3 y. (22) 



The second term of (12211 is irrelevant because it is zero on a surface of primary constraints, 
according to the condition (1201) . The rest of the PBs being known (see (fT6l) and (TT/Tj) ) allows 
one to just read off the functions a a (x, y) from (1221) and obtain an explicit expression for 

Gfo) (x) = -x 0<7 ~ ^oo,o^ 0,J + 9^,09^ - \g uk 9^ </>°° - (23) 

— 2 — -^=/ mn pfc7T — — + I m kpqQaBl — F^T^ 1 * 7 '' ^ I 0° 
0°° 0°° / 

/ Om \ 

mn,k ~r Qkm,n 9kn,m) I v 

This completes the calculation of the generator (T2T]) . 

V. TRANSFORMATION OF THE METRIC TENSOR 

Transformations of fields can be found by calculating the following PB 

Sg^ = {G, g^} = {kGfo + ^, o Oct , g^} ■ (24) 

Let us compare the result of this PB with ([1]) which we present in equivalent but more 
suitable form for comparison with (1241) 

^(diff)g^u = - + g al3 {g^p," + 9vp,n - g^,p) £<*■ (25) 

For the time-time component, g 0Q , (1241) gives 



figoo = {G,goo} = (Co-^fo) + C<r,o0 CT ) 

so that non-zero contributions come only from those terms proportional to the primary 
constraint (fP° 

Sgoo = ^oo,o^ 0<T + 9Gk,o9° k ~ ^g oh 9oo,kj ~ Co,o- (26) 
Putting fi = v = and partially separating the space and time indices in Eq. ff25l) . we have 

5(diff)goo = —2^o,o + g a0 goo,o£cT + g uk (2^,0 — goo,k) 

which is equivalent to (125]) up to a numerical factor 

2{^oo,C} = <5(di//)0oo (27) 

that can be incorporated into a parameter £ CT — > 2£ CT . 

For the transformation of the space-time components, g^k, we need the corresponding 
part of the generator depending on (p 0s 

1 1 1 / 1 g ap g am \ 

Sgok — ~2^0,k — 2^,0 — -jia (^—2 ^_ I mnpk TC mn — + Imk pq gal3,t-^-A (Pq) a j (28) 
+ o^c I <7 "POO,*: + 2g na g0n,k + g^ {.9mn,k + <?fcm,ra — gkn,m) 



In the first bracket of (1281) we substitute momenta in terms of the metric (j3J) so that using 
(JHj) we obtain 



qCct ^qq ( fl 1 gmk,n 9 gak,m 9 9am,k^) 



k 9 ~ 

and finally, after a simple rearrangement of terms in (|28|) . we have 



Sgok — — 2^o,fc — 2^ fc '° + 2^9° a 9m,k + 2^9™ (~9ok,n + <70n,ft + 9nk,o) ■ (29) 

This is equivalent to the transformations of the corresponding components of ( 1251) up to the 
same numerical factor that occurs in (127j) . 

For the transformation of the space-space components, g nm , we need to keep only terms 
in G with 7r pq dependence, so that 



5g nm = {G,g nm } = — ^ - ^—i^^—jl* 

The second term produces contributions proportional to the primary constraints and 
can be neglected on the surface of primary constraints. In x 0fJ we need only momentum- 
dependent terms appearing in (|T3|) . We obtain then 



$9nm 2^ n ' m 2^ m ' n ^ a ^a ^ e<? ^nm 2^™™ e ^ ^Pl,' 

>i / i \ 
I p £ I t Ott 11 ^ -\- T A{ c d)0nuk I 

g m 2 V V^g J 
After substitution of 7r ab from (j3J) and using flfJJ), and f|T4l) . we have 



$gnm 2^ n > m 2^ m i n 2^ ^ (.90m,n 90n,m 9nm,o) Q {.gpm,n gpn,m 5ntn,p)) 



(30) 



nq,m gmq,n gnm.q) 9 ^ (50m,n gOn.m ~ 9nm,o)) 



which is again equivalent with (1231) up to the same numerical factor of 2. 

The transformations of the components of g^ v in (1261) . fl29|) . fl30|) can be combined into 
one covariant expression (125]) or ([1]). This completes the proof of our original statement that 
the Hamiltonian formulation gives the same result for the gauge invariance of GR as the 
Lagrangian formulation [15( , as it should. Moreover, this also demonstrates that we obtained 
the consistent Hamiltonian formulation of GR because the transformation ([T]) derivable from 
it is generally covariant. 



VI. CONCLUSION 



We would like to point out some peculiarities of the Hamiltonian formulation of GR. 
First of all, the transformation (1317]) is valid on the surface of primary constraints. In 
ordinary gauge field theories (e.g., Maxwell or Yang-Mills) this does not happen, because 
the right-hand side of (I2"0"j) equals zero exactly. Another peculiarity of GR is the PB among 
primary and secondary constraints is not zero by (TT6j) . Such deviations from ordinary 
gauge theories actually can be expected. In ordinary field theories, equations of motion are 
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exactly invariant under a gauge transformation, whereas the Einstein equations of motion, 



G^u = Rfiu — \gnvR = 0, transform into a combination of the same equations [121 ]. i.e. are 
invariant only on-shell 

fi(diff)Gfj,v = —^ p Gfj, v>p — ^ffivp — ££G> P . (31) 

In particular, because of this, the "crucial test" of having off-shell closure of the constraint 
algebra [22J seems as unreasonably strong for GR. Moreover, it is impossible to expect 
off-shell closure for generators of a transformation which by itself produces only on-shell 
invariance of the equations of motion ( 13TT) . 

Going back to the "spiritual" statement of Hawking, we can make the conjecture or 
possibly draw the conclusion that if diffeomorphism is derivable through the Lagrangian 
or Hamiltonian approach, then the spirit of GR is "alive" , as the mathematical essence of 
the spirit of GR is to retain general covariance and that diffeomorphism can be derived 
from its structure. Vice versa, if diffeomorphism cannot be obtained in the Lagrangian or 
Hamiltonian formulation, then such a formulation is not the equivalent to GR. In particular, 
our Hamiltonian formulation of GR which does not resort to any change of variables or 



modifications gives the same diffeomorphism invariance as the Lagrangian approach [15 ]. 
On the contrary, the ADM Hamiltonian formulation gives different transformations which 
can be presented in the form of diffeomorphism only after field- dependent and non-covariant 
redefinition of parameters [lpj ]: 

_L _ (_ n 00\ -1/2.0 k _ p k _y_t0 /oo\ 

e ADM - \ 9 ) £ ADM — Uiff qo W/- ^> 

Moreover, the equivalence of the Lagrangian and Hamiltonian formulations implies that 
the ADM Lagrangian should give the same field-dependent redefinition of gauge parameters 



(152"]) if treated by the Lagrangian method of [151 ] . and it is the expected result [23j. The 



equivalence of results in both formalisms for asymptotically flat space-times was discussed 



in 



24]. 



In this paper we demonstrated that the consistent Hamiltonian formulation of GR can 
be obtained by considering the metric tensor as the canonical variable and not performing 
any change of variables. Following the Dirac conjecture and applying the Castellani proce- 
dure, we derived the gauge generators from the full set of first-class constraints. From these 
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generators, without redefinition of gauge parameters, we explicitly derived diffeomorphism 
invariance in the Hamiltonian formulation of GR for the first time. In the Hamiltonian for- 
mulation of GR, the covariance is not manifest but the final result on gauge transformation 
is presented in manifestly covariant form as in the Hamiltonian formulation of ordinary field 
theories. Our result is another illustration of the resiliency of Einstein's General Relativity 
and any approach, if correctly used, cannot violate its symmetry and leads to the generally 
covariant transformation ([1]). 
Note added in proof 

After completion of this work it came to our attention that, by a different from [3] method 
which was used in |lo| . the "specific metric-dependent diffeomorphism" f)32p of the ADM 



formulation was also obtained and extensively discussed in 



25|. 
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